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Abstract 

We study the effect of non-vanishing masses and mixings among neutrino 
flavours on the detection of neutrinos from stellar collapse by a water Cerenkov 
detector. We consider a frame- work in which there are four neutrino flavours, 
including a sterile neutrino, whose mass squared differences and mixings are 
constrained by the present understanding of solar, atmospheric and laboratory 
neutrino detection. We also include the effects of high density matter within 
the supernova core. Unlike in the three flavour scenario, we find that the num- 
ber of events due to the dominant process involving electron-antineutrinos 
changes dramatically for some allowed mixing parameters. Furthermore, con- 
tributions from charged-current scattering off oxygen nuclei in the detector 
can be considerably enhanced due to flavour mixing. We also present a com- 
parison between the two possible scenarios, namely, when only three active 
neutrino flavours are present and when they are accompanied by a fourth 
sterile neutrino. 
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I. INTRODUCTION 



In a recent paper |IJ] (hereafter referred to as I), we discussed in detail the signatures of 3 
flavours of neutrinos from stellar collapse. The analysis was confined to Type II supernovae 
(which occur when the initial mass of the star is between 8-20 solar masses). Based on 
the work of Kuo and Pantaleone 0, where they include mixing among all three neutrino 
flavours, we found that the mixing between neutrino flavours leaves non-trivial signatures in 
the detector. The main conclusion of the paper was that the effect of mixing is to produce a 
dramatic increase in the events involving oxygen targets ||. These will show up as a marked 
increase in the number of events in the backward direction with respect to the forward peaked 
events involving electrons as targets (more than 90% of which lie in a 10° forward cone with 
respect to the supernova direction for neutrinos with energies E v > 8 MeV). In the absence 
of any mixing, there will also be a few events in the backward direction due to CC scattering 
on oxygen targets. Furthermore, the observed T7 e p events are the largest in number as well 
as least sensitive to the mixing parameters within the three flavour scenario. Hence they 
provide a direct test of the supernova models. Since the angular distribution of these events 
is approximately isotropic [Kj, they may be used to set the overall normalisation. 



The above analysis was done assuming the standard mass heirarchy necessitated by the 
solar and atmospheric neutrino observations ||. In the meantime, several authors have 
looked at the possible signatures of neutrinos and antineutrinos from supernova collapse 
for Super Kamiokande and Sudbury Neutrino Observatory(SNO). Dighe and Smirnov || 
have looked at the problem of reconstruction of neutrino mass spectrum in a three flavour 
scenario. These authors as also Chiu and Kuo |7] have compared the signatures in the 
standard mass heirarchy and inverted-mass heirarchy. While these papers and I incorporate 
the constraints from solar and atmospheric neutrino observations, the important question 
of the mass limits that may be obtained from the observation of time delay has also been 
analysed by Choubey and Kar || (see also the review by Vogel |§). For a recent review 
which also discusses aspects of locating a supernova by its neutrinos in advance of optical 



observation, see ref. 10 



In this paper, we continue the main theme of paper I, and extend the analysis to include 
a fourth sterile neutrino as well as discuss some aspects of three flavour mixing not discussed 
in detail in I. The main motivation to extend the analysis comes from the fact that such a 
sterile neutrino is probably required in order to explain the recent results from the LSND 
collaboration [|ll|]. It is well known that an explanation of the solar and atmospheric neutrino 
puzzles involve two very different scales of mass-squared differences. As such, both of them 
cannot be explained unless one invokes at least three neutrino flavours. However, the scale 
required in order to understand the results from the LSND collaboration is different from 
both atmospheric and solar neutrino puzzles. It has therefore become necessary to introduce 
an extra- neutrino species, u s , which however must remain sterile due to the LEP constraint 
on the Z-width. 

The main feature of the analysis as emphasised in I — that the charged current (CC) 
events on oxygen nuclei (which show a preference to be in the backward direction) — is 
preserved both in three and four flavour scenarios. We also find that the signatures for three 
and four flavours are very distinct for some parameter ranges in the dominant isotropic 
events caused by electron antineutrino interactions with protons in water. 
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The analysis presented in this paper follows closely the earlier analysis presented in I 
and therefore we will only reproduce the main outline and refer to I for some details. In 
Sect. II we give an outline of the framework and the mixing matrix as also matter effects 
on mixing. We then use this to obtain expressions for the fluxes reaching the detector in 
Sect. III. Sect. IV contains an analytic discussion of the signatures in the case of extreme 
adiabatic and extreme non-adiabatic mixing. In Sect. V we write down expressions for the 
event rates due to interactions between different neutrino flavours and scattering targets in 
the detector. Numerical results are presented in Sect. VI both for the spectrum and the 
integrated number of events. We present a summary and list our conclusions in Sect. VII. 
Appendix A gives details of the choice of mixing matrix that we have used in the paper. 



II. MIXING IN THE PRESENCE OF HIGHLY DENSE MATTER 

We briefly discuss mixing among four flavours of neutrinos (or antineutrinos) and com- 
pute the neutrino survival and conversion probabilities. 

Unlike in the three flavour scenario, the addition of a sterile component forces us to 
prescribe a heirarchy of states in the mass eigenstate basis. If neutrino oscillation is the 
mechanism for the result from LSND, it is clear that the conversion of — > v e is governed 
by a mass scale in the range of 0.1 eV 2 to 1 eV 2 . As shown in Fig.l, we choose two doublets 
separated by this mass scale. In the lower doublet the mass-squared difference is given 
by the appropriate scale for the solution of the solar neutrino puzzle {5m% < 10 -5 eV 2 ). 
Analogously, in the upper doublet it is the atmospheric neutrino mass scale(5m\ TM ~ 
10~ 3 eV 2 ) that plays the crucial role. In principle the sterile neutrino may be in either of 
these doublets. However, if one believes that the atmospheric neutrino solution is driven by 
the conversion of — > u T , then the sterile neutrino should be in the lower doublet accounting 
for the solar neutrino deficit. While we take the sterile neutrino in the lower doublet, the 
final results are not crucially dependent on where we position the sterile neutrino as long as 
the two doublet scheme is adhered to. 

The mixing matrix, which relates the flavour and mass eigenstates in the four flavour 
scenario has six angles and the CP-violating phases which we do not consider here. As is the 
convention, we denote the mixing angle in the lower doublet by u and the mixing angle in the 
upper doublet by ip. The resulting complicated mixing matrix involving six mixing angles 



is greatly simplified by application of the chooz |T2| constraint. The chooz result implies 
that sin 9 < e, where 9 is either the (13) or (14) mixing angle and e < 0.16. We will therefore 
replace these angles by their maximum possible values allowed by chooz. Furthermore, we 
choose sin ^13 ~ sin 9\4 ~ sin #23 ~ sin #24 ~ £■ (The oscillation probabilities we are mainly 
interested in depend on the v e survival and conversion probabilities. These involve only #13 
and 9 14 mixing angles. As such #23 and #24 may be kept arbitrary and have been fixed to e 
for convenience. Note that the CHOOZ constraint corresponds to the replacement sin — > e 
in the three flavour case analysed in I, where is the (13) mixing angle to which the chooz 
constraint applies. This choice of hierarchy and mixings is consistent with already known 
data from various laboratory and atmospheric neutrino experiments. For more details on 
the constraints on mass squared differences and mixing angles from various experiments, see 
Appendix A. Then the flavour eigenstates are related to the four mass eigenstates in vacuum 
(for both neutrinos and antineutrinos) through a unitary transformation, 
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where the superscript v on the r.h.s. stands for vacuum. Within the two doublet scheme 
the 4x4 unitary matrix, U v , may be written as 
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where = sinu and c w = coscj, etc. The angles uo and ^ can take values between 
and 7r/2. The mixing matrix given above is unitary up to order e 2 . Since e is a small 
parameter, we do not need to go beyond this order. Later on we will show that the survival 
and oscillation probabilities that we need are such that they do not depend on the mixing 
in the upper doublet, namely the angle ip. For all practical purposes it can be set to any 
value, in particular, zero. 

The masses of the eigenstates in vacuum are taken to be /ii, /i2, At 3 and /i4. In the mass 
eigenbasis, the (mass) 2 matrix is diagonal: 
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= /i 2 1 + AM 2 , 

where the mass squared differences are given by $21 

Without loss of generality, we can take 8 2 i, 5 3 i and 5 4i to be greater than zero; this 
defines the standard hierarchy of masses consistent with the range of the mixing angles, 
as specified above. Neutrino oscillation amplitudes are independent of the first term so we 
drop it from further calculation. In the flavour basis, therefore, the relevant part of the mass 
squared matrix has the form, 
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U v AMg u^. 



(4) 



where U v is the mixing matrix in vacuum. 



A. Matter effects for neutrinos 



The relevant matter effects may be included by a modified mass squared matrix, 

AM 2 = AM 2 + AM a, (5) 
where the matter effects are included in, 
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with the Ai(r) given by 



and 



A x {r) = V2G F N e (r) x 2E 



A, 



V2 G F NJr) x E 



(7) 



(8) 



which are proportional to the electron number density, N e (r), and neutron number density 
N n (r) respectively. We set N e (r) ~ f e N n (r), where f e is the electron fraction which is 
less than one, in the supernova core. Here r is the radial distance from the centre of the 
star. Note that we have subtracted the neutral current contribution to the active flavours 
in the interaction term. Since the sterile neutrino does not interact at all, we have added 
and subtracted the neutral current contribution; hence the appearance of A 2 in the mass- 
squared matrix in matter. The detailed modifications due to matter effects are discussed in 
Appendix A of I. 

As noted in I, the maximum value of Ai occurs at the core and is approximately 2 x 
10 7 E eV 2 , where E is the neutrino energy in MeV. A 2 is somewhat smaller. The modification 
due to the matter dependence is similar to the case of solar neutrinos, although, unlike in 
the case of solar neutrinos, all active flavours are produced in the supernova core. 

It is clear that the mass squared matrix is no longer diagonalised by U v in the pres- 
ence of matter; we therefore diagonalise AM^ in order to determine the matter-corrected 
eigenstates. The value of Ai for energetic neutrinos (of a few MeV to tens of MeV) 
in the core is several orders of magnitude greater than these mass-squared differences. 
The eigenvalue problem may thus be solved perturbatively, with the following hierarchy: 
Aj(core) ^> <5 4 i > S 3 i ^> S 2 i- As a result, the electron neutrino undergoes three Mikheyev, 
Smirnov and Wolfenstein (MSW) resonances [ 13| . In Fig. 2 we show schematically the level 



crossing pattern for neutrinos in the presence of matter. Note that in this case there are 
several level crossings because of the presence of the sterile neutrino. 



1. The adiabatic case 

Following the method outlined in I, we compute the eigenstates at production (almost all 
the active neutrino flux is produced in the core). Unless e is extremely small, the propagation 
is adiabatic as may be shown by computing the Landau-Zener(LZ) probabilities at each 
resonance. While we include LZ probabilities in the actual numerical computation, we 
will first state the relevant probabilities in the case of adiabatic propagation of neutrinos 
produced in the supernova core. The average transition probability between two flavours a 
and (3 is denoted by P a p where a, (5 — e, /x, r, s. These are obtained by observing that almost 
all the v flux is produced in the highly dense core. Here Ai « 10 7 eV 2 , so one may take the 
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extreme limit in which the density is infinite. Then the electron neutrino is produced as the 
highest mass eigenstate, that is, 

We)"" = 

A similar conclusion holds for the sterile partner in the lower doublet. In the case of adiabatic 
propagation the mass eigenstates produced in matter remain the same in vacuum. Therefore 
to a good approximation, the probabilities (of survival or transition into each other) for these 
two flavours, v e and v^, are just the overlap of the highest two mass eigen states with the 
corresponding vacuum eigen states, and are small, as can be seen from the relevant entries 
in the matrix, U v , in eq. (0): 

P=P=P=P=e 2 (9) 

1 ee 1 es 1 se 1 ss c > 

independent of other mixing angles and the mass squared differences. Note that e 2 is con- 
strained by the limit set by the chooz collaboration and is small. Indeed in the purely 
adiabatic transition therefore the survival probability of the electron type neutrino is rather 
small. We also note that it is sufficient to know only these probabilities as long as the z/ M 
and v T are not separately detected. Furthermore, note that there is no initial flux for the 
sterile neutrino, but it can arise after oscillation due to mixing from other flavours. 



2. The non- adiabatic case 

A general non-adiabatic case is harder to discuss analytically. Since there are many 
crossings, one has now to consider non-adiabatic transitions at all these crossings. Because 
of the parametrization it is easy to see that non-adiabatic effects are introduced as a result 
of the values chosen for e and uo. The value of e determines whether non-adiabatic jumps are 
induced at the upper resonances while the value of u determines whether the non-adiabatic 
jump occurs at the lower resonance. This statement holds both for three and four flavours 
since in both cases the non-adiabaticity in the upper resonance(s) is controlled by e, apart 
from mass squared differences. 

Note that for a large range of e, allowed by the chooz constraint, the evolution of the 
electron neutrino is adiabatic. As a result the lower resonance does not come into the picture 
at all except when e is very close to zero, where the jump probability abruptly changes to 
one. The subsequent discussion is therefore relevant only when e is almost vanishing which 
is not ruled out by the known constraints except in LSND. The extreme non-adiabatic case 
occurs when the nature of transitions at the lowest resonance which involves the crossing of 
the first and the second mass eigenstates is nonadiabatic. Partial non-adiabaticity results 
when the transition in the lowest resonance is adiabatic while in the upper resonances it is 
non-adiabatic. The probability of jump at the lower resonance is given by Pl which is in 
general a function of to and the mass squared difference 5m\ 2 . In our calculations we have 
used the form discussed in the appendix of I (see also flU}]). 

The relevant survival and oscillation probabilities for electron-neutrinos in the non- 
adiabatic case are given by, 

P ee = (l-2e 2 )[(l-P L )sin 2 ^ + P L cos 2 ^], (10) 
P es = (l-2e 2 )[(l-P L )cos 2 ^ + P L sin 2 ^], (11) 
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P ss = P ee + 2e 2 (l-2P L )sin2cu, (12) 
P se = P es - 2e 2 (l - 2P L ) sin 2a;, (13) 

where u is as usual the vacuum mixing angle defined earlier. In the three flavour case the 
only relevant probability is, 

P ee = (l-e 2 )[(l-P L )sin 2 o; +P L cos 2 c4 (14) 

Note that, since e is small, the flux at the detector is entirely controlled by uo and Pl which 
is also a function of to. 



B. Matter effects for antineutrinos 

We now consider the case of V e propagation in highly dense matter. The only change in 
this case is that the matter dependent term in the relevant part of the mass squared matrix 
has the opposite sign (to that in Eq. (|7|)), that is, 

A(r) = -y/2 G F Ni(r) x2E;N l = N e , N 2 = N n /2. (15) 

The analysis goes through as in the case of v e propagation through matter. There are no 
Landau-Zener jumps to consider in this case since the resonance conditions are never satisfied 
unless the mass hierarchy is altered. The propagation is therefore adiabatic and the survival 
probability is obtained by simply projecting the eigenstate on to the flavour eigenstate 
in vacuum (at the detector). The antineutrino survival and transition probabilities for e and 
s flavour neutrinos are, 

P- e = (1 -2e 2 )cos 2 cj , (16) 
P- s = (l-2e 2 )sin 2 u; , (17) 
Pse = (l-2e 2 )sin 2 c<; + 2e 2 sm2c ( j , (18) 
P-= (1 -2e 2 )cos 2 cj-2e 2 sin2cj , (19) 

where u is as usual the vacuum mixing angle defined earlier. As before, we do not need to 
know the other probabilities for our analysis. 



III. FLUXES AT THE DETECTOR 

We briefly compute the neutrino flux at the detector in the presence of mixing and 
compare the three and four flavour scenarios. 

Following Kuo and Pantaleone 0, we denote the flux distribution, d0°/dP, of a neutrino 
(or antineutrino) of flavour a with energy E produced in the core of the supernova by P°. 
In particular we use the generic label P° for flavours other than u e and z/ e since 

F o = F o = F o = F o = F o (2Q) 

All these flavours are produced via the neutral-current (NC) pair production processes and 
therefore have the same flux for all practical purposes. However, the v e and V e fluxes are 
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different from each other and the rest since they are produced not only by pair production 
but also derive contribution from charged-current (CC) processes. Note that there is no 
production of sterile neutrinos at source. 

The flux reaching the detector from a supernova at a distance d from earth is reduced 
by an overall geometric factor of l/(47r<i 2 ). Apart from this, there is a further modification 
of the observed flux due to oscillations in the presence of matter. The flux on earth, in 
the various flavours, is given in terms of the flux of neutrinos produced in the core of the 
supernova by, 

F = P F° + P F° + P F° 

= P ee F° e + (1 - P ee - P es )F° x , (21) 

where we have made use of the constraint J2p Pa/3 = 1- Since and v T induced events cannot 
be separated in water Cerenkov detectors, their combined flux on earth may be written as 

= (P ee + P es + P se + P SS )F° + (1 - P ee - P se )F° , (22) 

and 

F s = (1 - P ss - P se )Fl + P se F Q e . (23) 

Note that a part of active neutrino flavours are converted into sterile due to mixing. The 
total flux is however conserved since, 

F e + 2F X + F S = F° e + 2F° . 

For comparison, the three-flavour case may be obtained by setting all the conversion 
probabilities involving the sterile neutrino to zero, and setting P ss = 1. 

IV. ANALYTICAL RESULTS FOR FLUXES 

Without any numerical estimation, or choice of supernova model, it is possible to analyse 
the adiabatic case as well as the extreme nonadiabatic case and obtain gross features of the 
mixing for both 3 and 4 flavour mixings. We shall analyse the mixings in the neutrinos and 
antineutrinos separately. We begin with the adiabatic case. 

A. Neutrino fluxes: Adiabatic case 

For adiabatic propagation, Landau Zener jumps must be negligible. This in turn means 
that the larger the mass difference, the smaller the mixing angle one can accomodate. We 
shall assume that this is satisfied for the case under consideration here. 

The results for the neutrino fluxes are given in Table |. It is seen that the fluxes are 
independent of the (12) mixing angle, u, in all cases. 
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1. 4~ flavours 

For the 4- flavour case, it is striking that the electron neutrino flux and the sterile neutrino 
flux (not observable in the detector) are equal and completely deplete the original spectrum 
in v^ T . The i>^ T flux at the detector is entirely made up by almost complete conversion 
of electron neutrinos. From the supernova models, it is reasonably well known that the 
spectrum of u e is colder than jT . In effect therefore the mixing interchanges the hot and 
the cold spectra in the supernova, with about half of the hotter spectrum lost in the sterile 
neutrinos. 

2. 3-flavours 

The 3-flavour v e flux is not very different from the 4-flavour one: about half of the 
hotter spectrum is lost and reappears as the electron neutrino spectrum in the detector 
(see Table |). There is no sterile neutrino; hence the v^ T spectra are not further depleted, 
in contrast to the 4-flavour case. However, this difference in v^^ spectra between 3 and 
4 flavour mixing is not easily observed since water Cerenkov detectors see very few events 
from flavours other than electron type neutrinos. The evidence for physical loss of flux into 
a sterile channel may be observed from neutral current events in detectors such as SNO. 

In both 3 and 4 flavour cases, there is mixing into the v e spectrum of the hotter u x 
spectrum. As we have seen in I, this admixture of hotter spectrum is signalled by a strong 
increase in the number of backwardly peaked events due to the opening up of the v O charged 
current channel in the detector, and should be a clean indicator of z/ e -z/ x mixing. 

B. Antineutrino fluxes 

The basic analysis is similar to the case of neutrino fluxes given in equations (|2l|) , (|22|) and 
(P3|). The changes occur only at the level of substituting for the appropriate probabilities, 
given in the previous section. These are different for neutrinos and antineutrinos. 

Unlike in the neutrino sector, predictions in the anti- neutrino sector require inputs on uo. 
Within the schemes considered in this paper, this is obtained only from an analysis of the 
solar neutrino deficit. Then there are two possible solutions for sin a; — the so-called small 
mixing angle (SMA) MSW solution, and also the large mixing angle (LMA) solutions. The 
LMA solutions arise both with matter effects included and in the so-called just-so solutions 
where there are no matter effects. At present, all three scenarios are consistent with the 
solar neutrino data although the mass-squared difference required in these cases is vastly 
different — 5m\ 2 ~ 1CT 5 eV 2 in the case of MSW solution with LMA, bm{ 2 « 1CT 7 eV 2 for 
LMA with vacuum (LMA-V) and bm\ 2 « KT 6 eV 2 for SMA solution. 

In any case, whatever the value of u, the propagation of antineutrinos always remains 
adiabatic (unless we choose an inverted mass hierarchy, which we do not consider here). 
The results for arbitrary uj are shown in Table ||. While there is hardly any admixture of 
the hotter T7^ T spectrum into T7 e (independent of u) in the 4-flavour case, the extent of this 
mixing depends on u in the 3-flavour case. 
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1. Small u 



It is easy to see from Table [TJ that for small u, the z/ e flux is almost unaltered, in both 
3 and 4-flavour mixing cases. These are then indistinguishable (and also indistinguishable 
from the no- mixing case). 



2. Large 10 

In 4-flavour mixing, the original V e flux is depleted by a factor of cos 2 u, while still 
receiving no contribution from z/^ iT (See Table [IT]). This depletion is seen at all energies 
and should be easily observable since the CC process V e p — > e + n has by far the largest 
cross-section in the detector; furthermore, it is approximately isotropic in distribution ||TU| . 



In contrast, a large value of to increases the contribution of v^ T in V e while proportionately 
decreasing the survival of the original V e in 3 flavours. Such a mixing will be observable at 
the higher energy end of the spectrum, where, however, the event rates are low. Hence, the 
behaviour of the isotropic events will be a clean signal of the number of flavours, provided, 
of course, that the original fluxes are accurately known. This depends on the models of 
stellar collapse. However, it may be possible to normalise the overall spectrum from the 
initial burst neutrinos (which are purely v e type and give a forward electron distribution in 
the Cerenkov detector). 

The V x flux is essentially unaltered, independent of u in the case of 4 flavours. With 3 
flavours, the original z/ x flux is depleted, and is compensated for by a proportional contribu- 
tion from the V e flux. Hence the results are reversed compared to the V e case. This signal, 
however, may be difficult to observe as there are very few events in this channel. 



C. Sterile neutrino fluxes 



We notice from Table | and Table [Tl] that in the four flavour case, there is actual loss of 
flux into the sterile channel. Independent of u, half of the hot v^ T spectrum is lost into the 
v s channel, while a portion proportional to sin 2 to of the v e flux is lost into the v s channel. 
From a detection point of view, this change, even if dramatically large, cannot be observed, 
unless through observation of neutral current events. 



D. Neutrino fluxes: Non-adiabatic case 

As we have seen, the various survival and transition probabilities in the non-adiabatic 
case depend on both e and u. For a large range of e the propagation is adiabatic. Hence 
non-adiabatic effects are important only for small values of e when the jump probability at 
the upper resonances abruptly changes to one. Then, the flux in this case is driven entirely 
by u, which also determines the size of the Landau-Zener jumps at the lower resonance 
through the value of the jump probability, P^. 
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1. Small u 



The neutrino fluxes at the detector, calculated using these probabilities, are summarised 
in Table III . For comparison we have also given the results for three flavours. It is clear from 
the table that if the mixing angle u is small, the v e flux at the detector is approximately 
given by, 

F e w P L F° e + (1 - P L )F° X , (3-flavours), (24) 
F e « P L F e ° (4-flavours) . (25) 

Hence the signals in the 3 and 4 flavour cases are drastically different: there is no mixing of 
the hot spectrum into the original v e spectrum in the 4 flavour case. 



2. Large to 

In the case of maximal mixing, u = 7r/4, we have 

F e «±[F° + F°] , (3-flavours), (26) 

F e « ^F e ° (4-flavours) . (27) 

Near maximal mixing, the LZ jump probability Pl plays no role at all. However, the v e flux 
at the detector again differs dramatically from three to four flavours (See Table [TTTJ ) . The 
three flavour result is also unlike the fully adiabatic case since we have an equal mixture of 
the cold and the hot spectrum due to mixing. Furthermore, the four flavour scenario differs 
from the four flavour adiabatic case quite dramatically. Unlike the adiabatic case, here there 
is no contribution from the hot spectrum, while the cold spectrum is depleted by about half. 

In short, the v e signal in the nonadiabatic case is very different not only between 3 and 
4 flavours, but also is very different from the adiabatic case. Here we have discussed the 
two extreme cases where the lower resonance undergoes or does not undergo Landau Zener 
jumps. We clarify the general case numerically below. 

In the case of anti-neutrinos, non-adiabatic transitions do not occur since there are no 
level crossings with the mass hierarchy assumed in this analysis. However, in the case of an 
inverted mass hierarchy that allows for non-adiabatic proagation in the antineutrino sector, 
the survival probabilities are identical to those given above. We will not discuss this situation 
further here. 

In the next section, we compute numerically the fluxes, for some preferred values of 
u and 5m\ 2 (obtained from an analysis of solar neutrino data) with the jump probability 
determined by these parameters. We will also predict event rates and the total number of 
events in the different cases we have discussed above. 



V. CROSS-SECTIONS AND EVENT RATES 

The basic quantity we are interested in is the distribution of events in the detector as a 
function of the energy of the detected particle. For a detailed discussion, for example, see I. 
Since we are interested in water Cerenkov detectors, the relevant interactions are, 
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Vi(pd + e" -> z^(z^) + e" , £ = e, /i, r ; 

I7 e + p -> e + + n ; (28) 

v e + 16 -f e" + W F , (29) 

77 e + 16 -> e+ + 16 N . (30) 

The oxygen cross-sections have been taken from Fig. 1 of Haxton || . All other cross-sections 
are well known. As the interactions on protons and oxygen nuclei are purely CC interactions, 
they involve only v e and V e . The interaction with electrons involves both CC and NC 
interactions for v e and V e and only NC interactions for all other flavours. The V e p cross- 
section is the largest, being proportional to the square of the antineutrino energy. In terms 
of total number of events, therefore, water Cerenkov detectors are mostly dominated by V e 
events. However, the different interactions in the detector have distinct angular distributions: 
The elastic neutrino-electron cross-sections are forward peaked, especially for neutrinos with 
energies > 10 MeV |1|JT|], while the proton cross-section is approximately isotropic in the 



lab frame. There is a slight excess of backward events for energies below 15 MeV, and slight 
excess in the forward direction at higher energies but the excess in either direction at 



the relevant energies are limited to few percent of the total number of events. Finally, the 
CC v e (v e ) cross-section on oxygen, although having a rather large threshold of 15.4 MeV 
(11.4 MeV) ||, increases rapidly with incoming neutrino energy and the number of events in 
the backward direction increases substantially with energy of the incoming neutrino above 
the threshold. 

The time integrated event rate, from neutrinos of flavour a and energy E, as a function 
of the recoil electron (or positron) energy, E e , is as usual given by, 

Here the index b refers to the time interval^. Note that the flux distribution, F a (b) includes 
the effects of mixing in the hot dense core and is a function of the time-dependent temper- 
ature, 1},. The index p refers to any of the various processes through which the neutrino a 
can interact with the target, t, in the detector. Here n t refers to the number of scattering 
targets (of either e, p or 16 0) that are available in the detector. The total number of events 
from a given flavour of neutrino in a given bin, k, of electron energy (which we choose to be 
of width 1 MeV) then is 

N ^ = i iE <it- (32) 

In the next section we give detailed numerical estimates both for the spectrum and the 
integrated events for some specific representative values of e and u, as constrained by other 
experiments. 



The constant factor At;, should also appear in the corresponding equation, Eq. (27), in I. 
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VI. RESULTS AND DISCUSSION 



As in I, we compute the time integrated event rate at a prototype 1 Kton water Cerenkov 
detector from neutrinos emitted by a supernova exploding 10 KPc away. Results for any 
other supernova explosion may be obtained by scaling the event rate by the appropriate 
distance to the supernova and the size of the detector, as shown in Eq. ([31]). We assume the 
efficiency and resolution of such a detector to be perfect. 

We use the luminosity and average energy distributions (as functions of time) as given in 
Totani et al. |18|], based on the numerical modelling of Mayle, Wilson and Schramm ||19 |. In 



a short time interval, At;,, the temperature can be set to a constant, T&. Then, the neutrino 
number flux is described, in this time interval, by a thermal Fermi Dirac distribution, 

for neutrinos of flavour a and energy E at a time t after the core bounce. Here b refers 
to the time-bin, t = t + bAt. Hereafter, we set the time of bounce, to = 0. The overall 
normalisation, Nq, is fixed by requiring that the total energy emitted per unit time equals the 
luminosity, C(a), in that time interval. The total emitted energy in all flavours of neutrinos 
is about 2.7 x 10 53 ergs, which is more or less equally distributed in all flavours. The number 
of neutrinos emitted in each flavour, however, is not the same since their average energies 
are different. 

It turns out that the results are not very sensitive to the time dependence of the temper- 
ature profile. For instance, using time-independent temperatures of 11, 16, and 25 MeV for 
the u e , V e and the v x thermal spectra (instead of the above time-dependent ones) changes 
the results by less than 5% in the antineutrino sector and by about 8% in the neutrino 
sector. However the difference between the average temperatures of the electron and other 
types of neutrinos is crucial to the analysis. 

In addition we impose the following known constraints on the mixing matrix in vacuum 
both for three and four flavour scenarios. Consistent with the chooz constraint, namely 
e < 0.16 (which we have already imposed at the level of the parametrisation itself), we 
set e = 0.087. Furthermore, the following constraints derived from solar and atmospheric 
neutrino observations are also taken into consideration: 

The constraint from the atmospheric neutrino analysis implies that the relevant angle 
7r/4), is near maximal and the relevant mass squared difference is of the order of 10~ 3 . 
Neither of these constraints directly enter our calculations except to determine whether the 
upper resonance is adiabatic or not depending on the value of e as constrained by the chooz 
findings. We consider both here. 

The more relevant constraint follows from the solar neutrino physics. Here there are 
three possible best fits to the combined data on solar neutrinos fl7| : 



1. sin 2 (2w) = 6.0 x 10- 3 ,5m 2 2 = 5.4 x 10~ 6 eV 2 (SMA). The small angle MSW solution. 

2. sin 2 (2w) = 0.76, 5m\ 2 = 1.8 x 10~ 5 eV 2 (LMA). The large angle MSW solution. 

3. sin 2 (2u;) = 0.96, 5m\ 2 = 7.9 x 10~ 8 eV 2 (LMA-V). The large angle vacuum solution. 
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While choosing only these values may appear restrictive, it will be seen that these values 
cover the typical ranges within which all the allowed changes take place. The numerical 
calculations are done by following the evolution of the mass eigenstates through all the 
resonances including the appropriate jump probabilities. 

We will now discuss the results for all these choices. We choose e = 0.087 for the adiabatic 
case and e ~ for the nonadiabatic one. Whether or not the Landau Zener jumps play a 
role is determined by the value of uj we choose. The interaction at the detector is mainly of 
three types: 

1. Isotropic events: These are by far the largest fraction of the events, and are due to 
V e p — > e + n. 

2. Forward peaked events: These are due to elastic scatterings of v e , V e , v x , V x on electron 
targets. 

3. Backward peaked events: These arise from CC scattering of u e and V e on oxygen nuclei 
in the target. 

Hence, angular information on the final state electron (e~ or e + ) will allow us to separate 
out the above three types of events. Note however that the extent of backward peaking is 
severely dependent on the spectral temperature . For the temperatures corresponding to 
the u x spectrum, we can assume the events to be totally in the backward direction. 

A. Total Number of events 

We give the time integrated number of events of the scattered electron, with energy, 
E e > 8 MeV (which is a typical threshhold for such detectors), in Tables 4 to 6, for the 
three choices of uj. These are the number of events for a supernova explosion at 10 Kpc for 
a 1 kTon water Cerenkov detector. 

The first column in all the tables corresponds to the base case of no mixing. The next 
two columns correspond to the adiabatic result (e = 0.087 for the three and four flavour 
case, while the last two columns correspond to the case where e ~ 0. The isotropic events 
(due to u e ) are insensitive to values of e chosen. The first table corresponds to small uj, 
where the lower resonance can involve Landau Zener jumps, while the next two correspond 
to large uj and are adiabatic at the lower resonance. 

It is expected that the integrated events are less sensitive to the details of models of 
stellar collapse, especially the energy distribution of the fluxes. Hence the predictions are 
more stable. 

1. e = 0.087 (Adiabatic case) 

The large backward excess due to mixing of a hot spectrum into the V e spectrum is a 
clear signal of mixing, irrespective of the number of flavours involved, as seen in Tables 4-6. 
Variations in the forward events may not be visible. However, the large depletion in the 
isotropic events with increasing uj (in the 4 flavour case) can distinguish between the 3 and 
4 flavour mixing, especially when uj is large (see Tables 5 and 6). 
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2. e ~ 



Here the lower resonance may be adiabatic or nonadiabatic, depending on the value of 
uo. Again here a large excess of backward events clearly indicates the presence of mixing, 
but this is seen only for the 3 flavour case (see Tables 4-6). The isotropic background 
from the V e p scattering is again sensitive to mixing and the number of mixing flavours, 
for large u (tables 5 and 6). Finally, unlike the previous case the depletion (enhancement) 
of forward events, independent of u, for 4(3) flavours may be observable, especially in the 
differential event rates, dN/dE e . It is therefore interesting to study the differential spectrum, 
which obviously has more information than the integrated events, to see whether this can 
be quantified further. We do this below. 

B. Event rates 

The event rates as a function of the scattered electron energy, E e , are displayed in Figs. 
3, 4, and 5. For detailed results in the 3-flavour case, see I. The solid lines in all the figures 
refer to the case when there is no mixing and serve as a reference. Both the three- and 
four-flavour results are displayed in each of these figures. 

In Fig. 3 we show the time integrated event rate per unit electron energy bin (of 1 MeV) 
as a function of the detected electron energy for the adiabatic case when e = 0.087. The 
dashed curves correspond to the case with mixing. Interactions with both electrons as well 
as oxygen nuclei occur. The event rates in both cases are shown. In both cases, mixing 
progressively enhances the high E e event rates. This shows up as an increase in both the 
forward and abckward peaked events. There is no perceptible difference between the 3- and 
4-flavour scenarios, independent of constraints on the angle u>. 

In Fig. 4 we show the electron neutrino spectrum for the case when e is small, in fact 
near zero. Hence the propagation near the upper resonances is nonadiabatic. Here, 3 and 4 
flavour mixing give drastically different results. For all u>, the forward and backward rates 
are depleted for 4 flavour mixing. However, these rates are small. The 3 flavour nonadiabatic 
mixing is indistinguishable from the adiabatic case (compare the left hand sides of Figs. 3 
and 4). 

Thus, an enhancement in the backward peaked events clearly indicates the presence 
of mixing. But there is no unambigious indicator of the number of flavours unless the 
parameters get restricted better through solar neutrino solutions. For this, we need to 
examine the electron anti-neutrino spectrum. 

In Fig. 5 we show the V e spectrum as a function of energy for three different choices 
of uj. (Here e = 0.087 but the results are essentially the same even if e is nearly zero.) 
The propagation is always adiabatic for the antineutrinos since there are no resonances to 
account for nonadiabatic jumps. The V e can interact with electrons, oxygen and protons in 
water and therefore contribute to forward, backward and isotropic events, the most dominant 
one being the approximately isotropic events due to interaction on the proton target (this 
cross-section is about two orders of magnitude larger than the others, for typical energies 
involved). For the SMA solution, there is no perceptible difference between the three and 
four flavour scenarios and the no mixing case, as already discussed. A perceptible difference 
occurs for LMA and LMA-V solutions. The intermediate energy (E e ~ 20 MeV) spectrum 
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is not depleted for 3 flavour mixing, while it is depleted for 4 flavour mixing (for the LMA 
solutions). Again, the 3 flavour isotropic spectrum is enhanced at larger energies (See 
Fig. 5); no such enhancement occurs in 4 flavours. Such an energy dependent pattern of 
event depletion and enhancement may help distinguish not only the 3 from the 4 flavour case, 
but also the mixing from the no-mixing solution. This pattern will survive, independent of 
uncertainties in the overall flux normalisations; however, it must be remembered that the 
event rate is very small at large energies (E e > 40 MeV). 

It is therefore seen that a combination of detection of both isotropic and forward events 
at low and intermediate energies (involving both T> e p and u e e scattering in water) can help 
distinguish the 3 and 4 flavour solutions for the case when e is close to zero. 

VII. CONCLUSIONS 

To summarise, neutrino mixing in the 3- and 4- flavour cases gives rise to very different 
event rates due to neutrinos (and antineutrinos) from stellar collapse interacting with a 
water Cerenkov detector. The neutrino sector is sensitive to the presence of mixing. This 
shows up typically as an enhancement of the higher energy events in the backward direction. 
This signal is not sensitive to the number of flavours. For a fairly large region of parameter 
space, the antineutrino events are sensitive to the number of flavours involved. Since this 
is the main signal of stellar collapse, the statistics should be good enough to make rather 
strong predictions from this channel. 

We would like to thank the organisers of the 6th Workshop on High Energy Physics 
Phenomenology (WHEPP-6), Chennai, January 2000, where a part of this work was done. 
We thank Sandhya Choubey, Srubabati Goswami, Kamales Kar, and Debashis Majumdar 
for discussions during this meeting. We would also like to thank Hari Dass for discussions. 

APPENDIX A 

We clarify here our choice of mixing matrix that we have used in the text. This arises 
from the requirement of consistency with various experiments. We discuss this, for the case 
of vacuum mixing, for both three and four flavour solutions. We shall set all CP violating 
phases to zero in what follows. 

For the sake of completeness, we briefly quote the results for the three flavour case, where 
all flavours are active neutrinos. We choose the mixing angles to be u, <fi and tp for the (12), 
(13) and (23) mixings. The electron neutrino survival probability is given by, 

P ee = (1 - 2c| S p _ 2c J4 s 2 + 2c 4 c 2 s 2 cos(2 .545 12 L/£) 

+244 [4 cos(2.545 13 L/£) + si cos(2.545 23 L/£)] , 1 " 1 

where L is measured in meters, S in eV 2 and E in MeV. 

We first assume a set of weak constraints, viz., that there exist two different scales for 
5: 612 < 10~ 5 eV 2 from solar neutrino data, and #13 ~ 10~ 3 eV 2 from atmospheric neutrino 
data. Hence <5i 3 ~ 5 2 3 = 5. 
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For the case of CHOOZ, we have L/E ~ 300 | 12| . Then the electron neutrino survival 
probability that can be applied in the CHOOZ experiment reduces to 



P ee = 1-4^1.27- > ( A2 ) 

independent of lu, where we have set sin0 = e. Hence the CHOOZ result constrains e to be 
small; in particular, e < 0.16, which justifies the approximation used in our parametrisation 



of the three-flavour mixing matrix. As a result one cannot accomodate the lsnd results JTT 
in a three-flavour frame-work. 

Hence it is necessary, for the sake of consistency of atmospheric neutrinos, CHOOZ and 
lsnd, to have at least one extra (sterile) flavour. 

We now discuss the four-flavour scenario within the mass hierarchy shown in Figure |l|. 
The same weak constraints as in the three-flavour case are assumed here. This results in a 
two-doublet structure for the four levels where the doublets are separated by the mass scale 
imposed by lsnd, namely, Si% ~ S14 ~ #23 — #24 ~ 0.1-1 eV 2 . In addition to these mass 
squared differences, there are now six mixing angles, 9 12 = 00, #34 = ip, 9\ 3 = 9i, 9u = 9 2l 
#23 = #3; #24 = #4 ; where ip is constrained to be maximal, ip = 7r/4, from atmospheric 
neutrino data. Hence the mixing matrix is 

U = U^U^UzU^ . (A3) 

We first discuss the constraint on P ee from CHOOZ in the 4-flavour case. Because L/E is 
large, of the order of 300, the oscillatory term involving Sis averages out to zero: 

= 0. (A4) 

The same holds for other mass squared differences of the same order of magnitude as #13. 
Furthermore, since S12 is small, we we can set cos(Si2L / (2E)) ~ 1. Therefore the survival 
probability of the electron neutrino is, 

P ee = 1 - - sin 2 29 2 - - cos 4 9 2 sin 2 20! + - sin 2 29 2 sin 2 29 x cos ( j , ( A5) 

2 2 2 \ 2E I 

where the scale £34 is set by the atmospheric neutrino data so that the oscillatory term is of 
order one. Note that the survival probability does not involve #3 and 64, as remarked in the 
text. 




Then the CHOOZ result, P ee = 1±0.04 [12] indicates that both 9 X and 9 2 should be small. 



We have therefore chosen these to be of the same order: sin 9i ~ sin 9 2 ~ e. Then 

P ee = l-4 e 2 + 2e 4 cos(^) , 
~ 1 — 4e 2 , 

independent of £34, L, and E. The CHOOZ null result therefore provides an upper bound on 
these mixing angles. 

Using all these constraints, the relevant transition probability for lsnd is given by 

P Me = 2e 2 [cos 9 3 - sin 9 4 sin 9 3 + cos 9 4 } 2 sin 2 ( j . (A7) 
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Here S is the separation between the two doublets. 

The lsnd result gives a non-zero value for P^ e , viz., P^ e = (2.6 ± 1.1) x 10 -3 [11|. This 
unambiguously fixes e to be different from zero even though the chooz result allows this 
value. Combining the limits from the two experiments, we have the constraint, 

clsnd < e < cchooz • {-A-o) 

The values we get in the four-flavour completely consistent set. We have used al- 

lowed values of the various parameters from the above in the analysis of supernova neutrinos 
in the text. 

Note that there are still no constraints on 63 and #4 which are the (23) and (24) mixing 
angles. Although we have set them to e in the text, the electron (and anti-electron) survival 
probabilities that we define are independent of these mixing angles (since they do not involve 
electron- type transitions). Hence this choice does not affect the results that we have obtained 
for these probabilities. Furthermore, we remark that the present generation of experiments 
is unlikely to constrain these angles significantly. A more detailed analysis of the allowed 
parameter space in the case of four-flavour mixing will be published elsewhere. 
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TABLES 



No. of flavours 


Neutrino flux at detector, Ff 


3 


F„ e = e'Fl + (1 - e 2 )F° 


4 


F, e =e 2 F° +(l-2e 2 )F° 


3 


2F X = (1 + e 2 )F° + (1 - e 2 )F° e 


4 


2F x = (4e 2 )F° + (l-2e 2 )F° e 


4 


F s = (1 - 2e 2 )F» + (6 2 )^ e 



TABLE I. Neutrino fluxes observed at the detector in the extreme adiabatic limit. 



No. of flavours 


Antineutrino flux at detector, Ff 


3 


F S = (1- 6 2 )c 2 F° + (4 + e 2 c 2 )F° 


4 


F,= (l-2e 2 )4F° + (2 e 2 )F° 


3 


2F, = (l + c 2 - e 2 4)F° + (4 + e 2 4)F, 


4 


2F* = (2 - 4e 2 )F° + +2e 2 (l - s 2u) )F B 


4 


F- s = 2e 2 F° + {sl-2e 2 {sl-s 2u} ))F- e 



TABLE II. Antineutrino fluxes at the detector. c w = cos a;, = sin a;. 



No. of flavours 


Neutrino flux at detector, Ff 


3 
4 


F e = (l-e 2 )[(l-P L )4 + P iC 2 ]F e ° 
+[ l_ ( l_ e 2 )((1 _p L)s 2 +PiC 2 )]F o 

F e = (1 - 26 2 )[(1 - P L )sl + F L c 2 ]F e ° + 2e 2 F° 


3 


2F X = [1 + (1 - e 2 )((l - P L )si + P L ci)]F» 

+ [ l_(l_ e 2 )((1 _p L)s 2 + p iC 2 )]jP 




4 


2F X = 2(1 - 2e 2 )F° + 2e 2 


1 + (1 - 2P l )s(2uj) 


F° 

- 1 e 


4 


F s = 2e 2 F° + {[(1 - P L )cl + P L 4] " 2e 2 [c 2 - 


-4-^LK + 2s 2 J]}F e 



TABLE III. Neutrino fluxes at the detector when non-adiabatic effects are introduced. Note 



that while the transition is assumed to be fully non-adiabatic at the upper resonances, it is con- 
trolled by the probability Pl at the lower resonance. 
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Events 


No Mixing 


Mixir 
3 flavour 


Lg(e = 0.087) 

4 flavour 


Mix 
3 flavour 


ing(e ~ 0) 

4 flavour 


Forward 


5.9 


7.6 


6.9 


7.1 


4.4 


Backward 


4.8 


27.3 


27.1 


21.5 


4.2 


Isotropic 


281.1 


282.0 


282.2 


281.3 


280.7 



TABLE IV. Total number of events (for two different choices of e) in the forward and backward 
direction, and isotropic events, respectively, for a scattered electron energy above 8 MeV. Here 
sin 2 (2w) = 6 x 10~ 3 , 5 2 m 12 = 5.4 x 10~6 eV 2 . 



Events 


No Mixing 


Mixir 
3 flavour 


Lg(e = 0.087) 

4 flavour 


Mix 
3 flavour 


:ing(e ~ 0) | 
4 flavour 


Forward 


5.9 


7.6 


6.6 


7.2 


3.8 


Backward 


4.8 


29.5 


26.2 


23.8 




Isotropic 


281.1 


306.7 


212.0 


306.1 


209.4 



TABLE V. Same as above for sin 2 (2u;) = 0.76, 5 2 m 12 = 1.8 x 10~5 eV 2 . 



Events 


No Mixing 


Mixir 
3 flavour 


Lg(e = 0.087) 

4 flavour 


Mix 
3 flavour 


ing(e ~ 0) | 
4 flavour 


Forward 


5.9 


7.7 


6.5 


7.0 


4.0 


Backward 


4.8 


30.7 


25.6 


21.8 


2.7 


Isotropic 


281.1 


320.8 


171.9 


320.3 


168.7 



TABLE VI. Same as above for sin 2 (2cj) = 0.96, 5 2 m 12 = 7.9 x 10~8 eV 2 . 
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FIGURES 



V 4 

8 m 5iTM 

V 3 



8 m LSND 





Sm ATM 


; 


5m| 



(a) 



(b) 



FIG. 1. The vacuum mass square differences in the 4 flavour schemes. Here v e and v s are 
predominantly mixed states of v\ and v<i while v„ and u T are that of v% and v\. The mixing 
between the lower and upper doublets has been chosen to be very small. Here S, ATM and LSND 
stand for the solar, atmospheric and LSND mass squared differences respectively. We have also 
shown the 3 flavour scheme for comparison. 




FIG. 2. Schematic drawing showing mass squares as functions of matter density in the 4 flavour 
scheme. Resonances occur at two different regions of matter density, the lower one at A « Sm 2 s 
and the upper one at A « ^w-lsnd- The upper one consists of 4 close resonances. The three flavour 
scheme is also shown for comparison. 
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FIG. 3. v e O and v e (for all flavours of v) event rates when the upper resonance is completely 
adiabatic. The solid lines represent the no mixing case and is plotted in all the graphs for compar- 
ison. The dashed lines are due to the effects of mixing. The oxygen events show dramatic increase 
due to mixing. Note that 3 and 4 flavour cases cannot be distinguished. 
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FIG. 4. v e O and ve (for all flavours of v) event rates when the upper resonance is completely 
non-adiabatic. The results depend upon the three possible solutions to the solar neutrino puzzle 
and are shown in the three panels, top, middle and bottom. The three flavour results are similar 
to the adiabatic case shown in Fig. 3 but the 4 flavour case shows suppression of the event rates 
in all cases. Here the different cases are distinguished by the extent of suppression. 
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3— flavour: Adiabatic 



4 — flavour: Adiabatic 
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FIG. 5. i7 e event rates for e = 0.087. The solid (dashed) lines are due to (no) mixing. While the 
3 flavour scheme shows enhancement of event rates at high energies due to mixing, the 4 flavour 
scheme shows suppression at lower energies. 
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